The Approach to the Thermodynamic Limit in Lattice QCD at 7^ 
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The expectation value of the complex phase factor of the fermion determinant is computed to 
leading order in the p-expansion of the chiral Lagrangian. The computation is valid for jj, < 
and determines the dependence of the sign problem on the volume and on the geometric shape of 
the volume. In the thermodynamic limit with Li ^ oo at fixed temperature I/Lq, the average 
phase factor vanishes. In the low temperature limit where Li/Lo is fixed as Li becomes large the 
average phase factor approaches one. The results for a finite volume compare well with lattice results 
obtained by AUton et al. After taking appropriate limits, we reproduce previously derived results 
for the e-regime and for 1-dimensional QCD. The distribution of the phase itself is also computed. 



I. INTRODUCTION 

Numerical lattice QCD at nonzero baryon chemical po- 
tential is obstructed by the sign problem: At nonzero 
chemical potential, /i, the phase factor of the fermion 
determinant 



2,g ^ det(g + At7o + rn) 
dct{D + ^1-/0 + m)* 



(1) 



invahdates a direct application of Monte Carlo methods. 
However, indirect methods have been devised to circum- 
vent the signproblem Q, i, i, i, i, @, 0, i, i, M, [HI, [12, 
[isl. [rn. Il6l|. Since these approaches only apply when 
the average of the phase factor is close to unity it is of 
considerable interest to understand when the fluctuations 
of the phase are mild and when they are severe. Since the 
measurement of the average phase factor on the lattice is 
plagued by the sign problem as well, it is imperative to 
understand the average phase factor analytically. 

In this paper we study the average phase factor ana- 
lytically within chiral perturbation theory. In particular, 
the approach to the thermodynamic limit will be ana- 
lyzed. 

Despite the absence of baryons, chiral perturbation 
theory has proved a vital tool in understanding lattice 
simulations at nonzero baryon chemical potential. The 
generating functional for the eigenvalue density of the 
QCD Dirac operator includes quarks with the opposite 
sign of the chemical potential [131 which therefore cou- 
ple to the third component of isospin. Usin g th is fact, 
the exact quenched [3 and unquenched [l9l. |20|| micro- 
scopic spectral density of the QCD Dirac operator were 
derived from a chiral Lagrangian. This result revealed 
that at nonzero chemical potential, the chiral condensate 
is related to the spectral density by a mechanism that is 
different from the Banks Casher relation [2l|: The dis- 
continuity in the chiral condensate is due to complex os- 
cillations on the microscopic scale [22[ in a macroscopic 
region of the complex eigenvalue plane. 

The microscopic limit is also known as the e-domain of 
QCD. Microscopic results for QCD can equally well be 



derived by means of chiral random matrix theory [23l . |2J] . 
This has the advantage that one may employ powerful 
random matrix methods such as orthogonal polynomials 
[25i . ,26. . .27. . i28| the replica trick [2^ or the sumpersym- 
metric method [s^l- For example, the unquenched mi- 
croscopic spectral density at nonzero chemical potential 
was first derived by means of random matrix theory [l^ , 
whereas the quenched spectral density at nonzero chem- 
ical potential was first obtained by means of the replica 
trick in combination with the Toda lattice equation . 

The recent computation of the average phase factor in 
the microscopic domain [U, shows that the average 
phase factor is suppressed exponentially with the volume 
when fi > m^/2. For such values of the chemical po- 
tential the quark mass is inside the cloud of eigenvalues 
of the Dirac operator and numerical lattice QCD simu- 
lations become exceedingly difficult. For smaller values 
of the chemical potential, the quark mass is outside the 
two dimensional domain of the eigenvalues, and the sign 
problem is less severe. 

In this paper we examine the character of the sign 
problem in the region < 171^^/2 and temperatures such 
that the use of chiral perturbation theory can be justi- 
fied. With /i < m7r/2 it was found in [sil, Is^l that the 
average phase factor remains nonzero in the microscopic 
limit where fiF-^^/V is held fixed as the volume is taken 
to infinity. For fiF-^ \/V ^ 1 the large volume asymptotic 
limit of the microscopic prediction is simply given by 



{^n^, = (1 



(2) 



(The quenched and the phase quenched average of the 
phase factor give identical predictions in this limit. Both 
are obtained by setting Nf = in the equation above.) 
This result suggest that unquenched lattice simulations 
in this domain are feasible. 

Here we examine whether the average phase factor re- 
mains nonzero for /i < to^/2 when we relax the micro- 
scopic constraints and approach the thermodynamic limit 
at fixed chemical potential. In order to do so we compute 
the average phase factor using the p-expansion of chiral 
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perturbation theory where 

p-l/L, m^^l/L, i-L^l/L, T-l/L, (3) 

and work to leading (one-loop) order. The previous 
calculation of the average phase factor was worked out 
[11] in the microscopic domain where m'^F^ ^ 1/V 
and n'^F^ ~ 1/V as the volume is taken to infinity. The 
new one-loop computation presented here includes the 
effect generated by the nonzero momentum modes of the 
Goldstone bosons. We keep explicitly the dependence 
on volume V = L^Lq and the ratios Lq/Li in order to 
study the approach to the thermodynamic limit. The 
new result bridges the gap between the microscopic pre- 
diction [3l|, [3^ and the parameter range typically used in 
lattice gauge theories. This allows us to compare the one- 
loop result for the average phase factor to lattice results 
by AUton et ai. Below the pseudo-critical temperature 
for chiral symmetry restoration the lattice results are in 
remarkably good agreement with the analytical predic- 
tions. 

The distribution of the phase itself (rather than the 
phase factor) also follows from the one-loop computation. 
We give the explicit form of the distribution of the phase 
to one-loop order in chiral perturbation theory. 

The paper is organized as follows. In the next section 
we present the general setup for computing the one-loop 
result for the average phase factor within chiral pertur- 
bation theory. The explicit one-loop result is derived in 
section IIIII This expression is evaluated numerically in 
section IIVI and the comparison to lattice data is made in 
section |Vl The effect of a finite box on the average phase 
factor is further discussed in section IVTl Section |VlI] con- 
tains the discussion of the distribution of the phase. We 
end with concluding remarks in section rVIIII 



II. THE AVERAGE PHASE FACTOR TO 
LEADING ORDER 

The average phase factor in the full theory is the ra- 
tio of two partition functions. A partition function with 
an extra fermionic quark as well as a conjugate bosonic 
quark divided by the usual QCD partition function 

{e'"')N, = (4) 
{dct^f+\D + M7o + m)/ detjP - njp + m)) 
{dct{D + n-fo + m)^f) 

Here we used that conjugate quarks correspond to ordi- 
nary quarks with the opposite sign of the chemical po- 
tential . With the usual setup of leading order chiral 
perturbation theory, see for example (s^ . l35j | . the free en- 
ergy is a sum of contributions from each of the Goldstone 
bosons. The contributions to the numerator of pions with 
no isospin charge cancel against contributions from the 
denominator. This leaves us with 

-Go(/x)) 



2i0 



where each Gq includes the contribution of two oppositely 
charged Goldstone modes. In order to get the combina- 
torics right, notice that the inverse determinant in (j4]) 



represents a conjugate bosonic quark. The charged Gold- 
stone bosons contain this bosonic quark in addition to 
one of the Nf + 1 fermionic quarks and are thus fermionic 
in nature resulting in an additional minus sign from the 
fermionic loop. 

Because of the sign problem one often studies averages 
in the phase quenched theory where the phase of the 
fcrmion is ignored 



{\det{D + ^ijo + m)\'^ 



(6) 



The average phase factor in the phase quenched theory 
is defined by 



|det(L»-HAi7o-Hm)|2) Zr+x- 



(7) 



To leading order in chiral perturbation theory this ratio 
is given by 



(8) 



Notice that this result coincides with the one-loop re- 
sult for the quenched theory obtained from ([5]) by setting 
N f = as well as with the result for partially quenched 
computations with dynamical quarks at zero chemical 
potential. 

In summary, to determine the average phase factor to 
one-loop order all that is required is the difference be- 
tween Go{fi) and Go{^^ = 0). We will derive this differ- 
ence in section Hm 



III. EVALUATION OF Go(/i) - Go(m = 0) 

As we are particularly interested in the approach to the 
thermodynamic limit we will consider a finite system with 
volume V = L^Lq. Hasenfratz and Leutwyler [34| worked 
out Go for /i = 0. This calculation was generalized to 
nonzero chemical potential in [33 |. For completeness, we 
repeat the main steps of the computation below. 

Consider a single charged Goldstone boson with charge 
2 in a box V = L^Lq. The one- loop contribution to the 
partition function is given by 



= Go(p)/2 



where 



exph^ ^log(pI + ml + {pko- 2z/i)2)], (9) 



Pkc 



2'Kka 

Pka = — p — , ka mteger. 



La 



(10) 



While Go (a*) is divergent the difference between Go(^) 
and Go(/i) for y = cxd is finite, that is 



Go(m) = Go(Ai)|y=oo + go{^J-) 



(11) 



with QoifJ-) finite. The po'integration contour in the sec- 
ond term of the difference 
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can be shifted by 2in if there are no obstructions from singularities. This is the case if 2^ < m^, so that the /j,- 
dependence resides entirely in go^n). The difference of the free energies in ([Sj and ([8]) is thus given by the difference 
of the finite parts 

Go(/i) - Go(m = 0) = g^{p) - .9o(a^ = 0). (12) 

This also shows that the average phase factor does not depend on the ultra-violet cutoff. The infrared nature of the 
average phase factor has been verified on the lattice jsH]. 

After several manipulations including Poisson resummation and Jacobi's imaginary transformation (the steps are 
given in detail in (3^ ) we find two equivalent representations of go{n) both valid for /i < m^/2, 

^e-^"^^V4.(-Q^e-Wo^-e-tl#-l). (13) 



and 



a=0 la 

where Iq. runs over all integers and we have introduced the length 

L^{L,Lif/\ (15) 

The first representation has the advantage that the chemical potential explicitly appears in the combination iiLq. 
(Note that Lq also appears in other places.). The second representation can be used to expand the result in a power 
series in and /x^. It also may be preferable for numerical evaluations. 



The one- loop result is valid for LAqcd ^ 1, mVY, ^ 1 
and both <C Aqcd and fi <C Aqcd- Finally, the 
condition /i < to^/2 has to be satisfied. The reason is 
that for /i > ?Ti^/2, the Goldstone fields have to be ex- 
panded about a rotated ground state. For m^L^ <C 1 
and y?L^ the dominant contribution is given by the 
zero-momentum term and the one-loop result reduces to 
(l2|). A small tti-kL and small ^.L expansion about this 
result was given in [s^ . The general one- loop result is 
also valid when m^L^ and /x^L^ are of order or larger 
than 1. In the next section we study the one- loop result 
numerically. 



IV. NUMERICAL EVALUATION OF THE 
ONE-LOOP RESULT 

In this section we study the one-loop expression for the 
average phase factor as a function of dimensionless com- 
binations of the parameters. The average phase factor 



depends only on three dimensionless combinations: 

rriTrL, fiL and — . (16) 
Li 

As the one-loop expression is only valid for /i < 
we will express our results as a function of the ratio 
/i/(TO7r/2). We will consider 3 different cases: 1) The 
thermodynamic hmit at fixed temperature where Lim^^ 
becomes large while Lo^tt remains fixed, 2) The ther- 
modynamic limit at low temperatures where Lim^^ grows 
for fixed asymmetry Li/ Lq^ and 3) The low temperature 
limit in a finite box where L^m-j^ increases for fixed Lim^^. 

Thermodynamic limit at fixed temperature: The 

thermodynamic limit at fixed temperature is obtained 
by letting the spatial extent of the box in units of the 
inverse pion mass go to infinity while keeping the tempo- 
ral extent also in units of l/?n^ fixed. In this limit the 
sign problem is acute for any nonzero value of the baryon 
chemical potential, as indicated by a vanishing average 
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FIG. 1: The average phase factor goes to zero as the spatial extend of the box (in units of 1/m^) becomes large as compared 
to the temporal extend. The approach to the thermodynamic limit is illustrated by varying 771.^ Li at fixed temperature l/Lo- 
Left: The (phase) quenched average phase factor. Right: The average phase factor with two dynamical flavors with the same 
mass. Top: m^Lo = 2. Bottom: m,r^o = 4. Note that the average phase factor starts dropping to zero when Li exceeds Lq; 
the spatial length Li needs to be larger than Lo for the effect of the pion loop to be present in the average phase factor. 




FIG. 2; The average phase factor when approaching the thermodynamic limit for fixed Li/Lo. Left: The (phase) quenched 
case. Right: Two dynamical flavors. Top: The temporal extent of the box is taken equal to its spatial extent, hence the 
j/-axis is also the inverse temperature in units of l/m,r. Bottom: The asymmetry is now set to Li/Lo — 2. Hence the I/m,r 
is also y/SLom^. In all cases, in the thermodynamic limit the average phase factor approaches a step function which becomes 
zero beyond 2/i — m,r. 



phase factor. In figure [T] we show the approach of the average phase factor to zero as Li increases. The two left 
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plots show the phase-qucnchcd prediction (which, as dis- 
cussed above, is equivalent to the quenched prediction) 
while the two right plots show the case with two dynam- 
ical flavors. Note that for Li < Lq the average phase 
factor is dominated by the static pion modes cf. ([2]) [5^ . 

Thermodynamic limit v^fith temperature going to 
zero: As is common practice in lattice QCD, we now fix 
the asymmetry of the box, Li/Lo, and vary the size of 
the box. A large box now automatically implies a low 
temperature. In this thermodynamic limit the average 
phase factor approaches unity as long as the chemical 
potential is less than half of the pion mass. For larger 
values of the chemical potential the average phase factor 
is zero. In figure [2] wc show the approach to this step 
function as a function of Li. The two top panels are for a 
cubic box (quenched and unquenched) and the two lower 
plots are for Li = 2Lq. The approach to the step function 
is slower when the asymmetry is larger than unity. 

Zero temperature in a finite box: Here we consider 
the limit where the spatial size of the box is fixed in 
units of I/ttItt as the temperature (in units of m,r) is low- 
ered. In the zero temperature limit where the asymmetry 
Li/Lf) ^ 1 the average phase factor again approaches the 
step function 6(771^^ — '2fJ.). This is not surprising given 
that this behavior was also found in the exact solution of 
one-dimensional QCD [s^l- 

In figure [3] we show how the average phase factor 
(quenched, top, unquenched, bottom) behaves as a func- 
tion of and T in a finite box. 




FIG. 3: The average phase factor for rriTyLi =3. It ap- 
proaches zero with increasing temperature. Top: The (phase) 
quenched case. Bottom: The Nf = 2 case. Notice that the 
zero temperature Umit is the step function ^(m^ — 2fi). In 
this limit the time direction is much longer than the spatial 
ones, and as expected, the step function found is consistent 
with the result from one-dimensional QCD [37| . 



V. COMPARISON TO LATTICE DATA BY 
ALLTON ET AL.. 

The main limitations of lattice studies of QCD at 
nonzero chemical potential are the phase fluctuations of 
the fermion determinant and it is natural to analyze them 
quantitatively 0, M, M, M, M, SH, \M, S, H H H 
[43. In this section we compare the one- loop results ob- 
tained above to the lattice data of AUton et al. [l^l ■ 

In [l^l the response of the QCD partition function to 
a baryon and isospin chemical potential was measured 
at zero value of the chemical potential. The response to 
second order in the baryon and isospin chemical potential 
was given in terms of two numbers C2 and C2 respectively 
(see table 3.2 of They are related to the quark and 

isospin susceptibility at = according to 



C2 



Xq 

2T2' 



XI 
2r2' 



(17) 



To second order in the chemical potential, the mea- 
sured average phase factor (in the phase quenched the- 
ory) is given by 



\e ;iat 



(18) 



,(c2-c|)(L./Lo)'(2M/m^)'('n./Tc)V(T/Te)V4 



where, in the second line, wc have expressed the result 
in terms of accessible dimensionless ratios. Note that 
the strength of the sign problem to lowest order in the 
Taylor expansion only depends on the coefficient of the 
off-diagonal susceptibility d^'^ = (c2 — C2)/4. 

The one-loop chiral perturbation theory result is ob- 
tained by matching the dimensions of the box and the 
chemical potential in units of the inverse pion mass to 
those of lip. To fix the scale we use the value m^/Tc = 
3.58 from [4^|. For temperatures below the critical tem- 
perature, Tc, the agreement is very good, see figure H] 
Since chiral perturbation theory is not applicable in the 
chirally restored phase the disagreement for T > Tc is as 
expected. Unfortunately, the isospin susceptibility was 
not calculated beyond second order in /i/ in (isj so that 
we cannot extract the average phase factor to higher or- 
der. The fourth order term is of particular interest be- 
cause the analytical one- loop result, displayed in figure 
m is well approximated by exp(— a/i^ — b^'^) where a and 
6 are positive constants. 

Beyond the critical temperature the severity of the sign 
problem decreases significantly. In the high temperature 
limit the difference from one of the average phase factor 
comes from terms of order (7^ and higher as can be in- 
ferred from [4^, [lOl- Indeed, the terms up to order ff^ 
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FIG. 4: The phase quenched average phase factor obtained 
using the data of (solid curves) compared with the re- 
sult from chiral perturbation theory (dashed curves). The 
temperatures are T/T^ = 0.76, 0.90, 1.00, 1.11 where Tc is the 
critical temperature at = 0. The agreement is good even 
close to Tc. Above the critical temperature the response to 
the baryon and isospin chemical potential becomes alike as 
bound states of quarks have melted. The prediction of chiral 
perturbation theory of course fails in this region. 



modes is ± 2/i (recall that these Goldstone modes are 
made out of a quark and a conjugate quark). As /i ap- 
proaches TO7r/2 the lightest mode becomes massless which 
invahdates the replacement the sum over momenta by an 
integral. 

VII. THE DISTRIBUTION OF THE PHASE 



are functions of /^^ so that the phase quenched and 
the unquenched partition function are identical up to this 
order. In lattice simulations by AUton et al. the quark 
and isospin susceptibility are very close for T > Tc- The 
physical implication is that bound states of light quarks 
are absent beyond Tc. 



VI. FINITE VERSUS INFINITE BOX 

In this section we compare the result for a finite 
box with its thermodynamic limit which has been used 
in the resonance gas model. 

In the resonance gas model it is usually assumed that 
Ljm,r 3> 1 such that the spatial momenta in the expres- 
sion for go can be integrated over instead of a discrete 
summation. This leads to the standard expression 



50 (/^) 



VmlT^ A i^2(%^ 



) , ,2^n 
-cosh( — ). 



(19) 



In addition, the resonance gas model includes heavier res- 
onances. Among others, it was applied [sH to the quark 
and isospin susceptibilities for temperatures beyond Tc- 
In Figure E] we compare the average phase factor in a 
box with finite spatial length to the result obtained using 
(fTO)) . We observe that the validity of the standard expres- 
sion (|19p depends on the chemical potential. This is not 
surprising. After all the mass of the charged Goldstone 



In addition to studying the average phase factor, it is 
natural to also analyze the distribution function of the 
phase [H, 112 itself. It is defined by 

P^^{9) = {S{9-9'))n^ (20) 

_ / dA\ det{D + /i7o + m)\^f e'^j'^' S{0 - 6'')e"^™ 
~ / dA\ det(D -I- njo + m)|^/e*^/«'e-'S'™ " 

The unquenched ^-distribution can be written 



p^,(0) = e^^/«p^^,/2+JV,/2.W■ 



^Nf/2+Nf/2' 



(21) 



where the phase quenched (^-distribution is defined as the 
average 

PNf/2+Nf/2' {6) = {S{0 - 0'))Nj,/2+Nf/2' (22) 

with respect to the phase quenched partition function 

ZNf/2+Nf/2' = (I dct{D + fijo + ™)rO- (23) 

This rewriting shows that the unquenched ^-distribution 
is complex for any nonzero value of the chemical po- 
tential. The complex nature of the unquenched 0- 
distribution resides entirely in exp{iNf9) - the other fac- 
tors in (pij) are real, positive and even. Despite its simple 
form, the effect of the phase is drastic. Integrating over 
8 we find that 



dee'^f'pN,/2+N,/2'{0) 



z 



Nf/2+Nf/2' 



(24) 
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which becomes exponentially small (~ exp(— V)) in the 
thermodynamic limit at fixed temperature. 

In principle, one can extract the unquenched partition 
function from the phase quenched one and the phase 
quenched ^-distribution using (|24p . Numerically it is, 
however, very difficult to handle the detailed cancella- 
tions. It only works if the width of ^-distribution is com- 
parable to 2'K/Nf so that the complex oscillations have 
little effect. Below we will compute the ^-distribution in 
a finite box using chiral perturbation theory and show 
that the width of the ^-distribution is rather of order 



becomes a simple Gaussian: 



A. ^-distribution to one-loop in chiral perturbation 
theory 

The distribution function of the phase, pjv/ (^), can be 
extracted from the moments of the phase factor. Just 
like the average phase factor the higher moments follow 
from the one-loop computation 



IN, 



- „n(Af/+n)(Go(M=0)-Go(M)) 



(25) 



This holds for both both positive and negative integers 
n. Using the replica trick [1^, H^, one can analytically 
continue n to noninteger values. For instance, 



INf 



,i(JV^ + i)(Go(M=0)-Go(M)) 



(26) 



We expect that the analytic continuation of (exp(2m6')) 
in n is valid when the quark mass is outside the support 
of the eigenvalues - as is always the case in this paper. 

In we evaluated the even Fourier components. If 
we assume that the odd Fourrier coefheients are given 
by the same expression we can simply express the delta- 
function in the definition of pMj{d) into a sum over the 
moments. Introducing the shorthand 



AGo = Go(/x)-Go(// = 0) 



(27) 



wc find 



oo 
n— — oo 

= lgiAf/e+iw|AGo g-me-ri^AGo/4 

TT ^ 

n— — OQ 

= ie«^/e+i^/^AGo^3(0/(2^), e-'^^"/^). (28) 

TT 

By a Poisson resummation this can be rewritten as 



1 



p^AO) = - ,.jV,g + ljV|AGoy,-(9+2n.)VAGo 



(29) 



valid for a compact phase angle 9 € [— tt, tt]. For a contin- 
uous phase angle, 6 G [— oo, oo] the distribution function 



PNf {0) 



1 



^{Nf/2fAGo„iNje-e''/AGo 



(30) 



The quenched as well as the phase quenched average are 
given by pNf=o{())- Notice that result ([50)) is consistent 
with the general form given in (|2ip . Since AGg is exten- 
sive, the width of the ^-distribution is of order \/V while 
its amplitude increases exponentially with the volume. 
Along with the fact that the distribution is normalized 
to one, this illustrates just how intricate the cancellations 
will be, and therefore how tough the sign problem will be 
to handle numerically. 

With ptvj (^) at hand it is straightforward to compute 
also the variance of the phase 



(31) 



Note that the result is independent of iV/ even though 
{(^'^)nj and (6')7Vf both depend on the numbers of fla- 
vors. This suggest that the variance of the phase can be 
obtained from the quenched theory. 

The Gaussian form of the phase quenched 9- 
distribution is in agreement with the numerical result of 
Ejiri [s^. Notice, however, that an error of order 
in the numerical determination of the width of the Gaus- 
sian will lead to an error of order one in the unquenched 
partition function, c.f. (j24p . Small non-Gaussian cor- 
rections can have a similar dramatic effect. We expect 
that higher order terms in chiral perturbation theory as 
well as effects from baryons will result in a non-Gaussian 
form. 



VIII. CONCLUSIONS 

The average phase factor of the fermion determinant 
has been computed and examined to one-loop order in 
chiral perturbation theory for quark chemical potential 
less than half of the pion mass. In the ordinary ther- 
modynamic limit at fixed temperature the average phase 
factor is zero for any nonzero value of the chemical poten- 
tial. If the temperature is taken to zero at fixed aspect 
ratio of the box, the phase factor remains unity when the 
chemical potential is less than m7r/2. This indicates that 
QCD at zero temperature has a mild sign problem for 
p < in the thermodynamic limit. It would be of 

interest to study this region on the lattice. In particular 
to examine if there is a spinodal line in this region. 

The one-loop prediction for the average phase factor 
is in agreement with lattice data below the critical tem- 
perature. Below Tc the one-loop prediction thus gives a 
direct way to estimate the strength of the sign problem 
for a given lattice volume V = L^Lq and quark mass. 

The distribution of the phase itself was also derived 
from chiral perturbation theory. Its simple Gaussian 
form is consistent with recent lattice simulations 1521. 
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The critical isospin chemical potential beyond which 
pions Bose-Einstein condense is expected to depend on 
the temperature. The effect of this shift on the strength 
of the sign problem is not included in the present paper. 
In lattice simulations such effect has been observed [SJI 
and it would of great interest to extend the present work 
in this direction. 
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